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Some new results on Duffie-type OTC markets 

BY ALAIN BELANGER*, GASTON GIROUX* AND NDOUNE NDOUNE* 

Universite de Sherbrooke 


Abstract: The extended Wild sums considered in this article generalize the classi¬ 
cal Wild sums of statistical physics. We first show how to obtain explicit solutions for 
the evolution equation of a large system where the interactions are given by a single, 
but general, interacting kernel which involves m components, for a fixed m > 2.We 
then show how to retain the explicit formulas for the case of OTC market models 
where the dynamics is more directly described by two (or more) kernels. 


1 Introduction 

After the publication of M. Kac’s work (1956) [8], there was a renewed interest 
for the results of E. Wild (1951) [15]. This interest was mainly focused on the 
random matching of a large population of particles forming a diluted Maxwell 
gas. Here we develop an approach inspired by this body of work. To do so, 
we start with a sequence of dynamical sets of interacting components, one for 
each integer N. For these dynamical systems we can show that when N is large 
the probability is very small that a component has interacted more than once, 
directly or indirectly, up to time t, with any other component. Thanks to this 
fundamental property, we can link the microscopic and macroscopic levels using 
results from the theory of continuous-time Markov chains. 

The Wild sum is a series construction which gives the solution of a given 
evolution equation in the statistical physics of gases as first appeared in the work 
of E. Wild [15]. Note that the classical expression of a Wild sum is described by 
binary trees. Inspired by these ideas, S. Tanaka [12] and H. Tanaka [11] defined 
an extension of Wild’s sum for solving certain non-linear differential equations 
of spaces of measures, so the expression of this sum is described by appropriate 
trees. However, the problem of showing the existence of these sums remains 
wide open in general. 

The recursive time relaxed Monte Carlo methods of Trazzi, Pareschi and 
Wennberg [14] are based on generalized Wild sums. However, the lack of explicit 
formulas for these sums constitutes a handicap for the efficiency of the above 
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methods as well as others also based on extended Wild sums (see [13], for 
instance). 

Carlen et al [4] obtain Wild sum formulas which are quite explicit for the 
solution of the Kac equation. Their binary trees are obtained, in the spirit of 
McKean, from commutator formulas for Lie algebras, leading them to group¬ 
ings of interaction trees. Consequently, our more general interaction trees are 
different form theirs even in the binary case. 

The aim of this paper is to propose a combinatorial formula for extended 
Wild sums which are solutions of certain evolution equations and more precisely 
in the context of interactions involving in components, to > 2. 

In section 3 of Belanger-Giroux [1], the explicit formulas for the Wild sums 
were used to obtain the convergence of the solution of the evolution equation 
to a steady state. This is one of the important applications permitted by the 
tractability of our explicit formulas. 

The article is organized as follows. In section 2, we introduce the types 
of combinatorial trees which are going to be useful in the expression of the 
solution of the evolution equation in terms of interaction trees. In section 3 
and 4 we consider interactions involving m components, for to > 2, and we 
suppose that the intensities of these dynamics have an adequate dependence 
on N. Our techniques enable us to obtain an explicit formula for the solution 
of the associated system of differential equations. In section 4, we show how 
to retain the explicit formulation of the solutions in the case of OTC market 
models described by two kernels. 

2 Combinatorial trees 


We assume that the reader is familiar with the basic definitions of trees. A 
rooted tree is a tree with a designated node called the root. A rooted tree in 
which the rooted node has one child is a planted tree. An to- ary tree is a rooted 
tree where each of its node is either a leaf (that is, it has no child) or it has 
exactly to children. The leafs are called external nodes and those nodes with m 
children, internal nodes. Note that we do not consider the root of the tree as 
an internal node. 

An ordered tree is a rooted tree in which the children of each node are assigned 
a fixed ordering. 

A rooted tree is called an (m, l)-ary tree if each internal node has either one 
child or exactly m children. In this article, we will work with ordered to- ary 
trees and ordered (to, l)-ary trees. 

Let A„ denote the set of m-ary ordered trees with n internal nodes. Each 
tree in A„ has (m — 1 )n + 1 leaves and each tree can be obtained by adding an 
internal node on a leaf of a tree in A„_i (taking into account the order). Hence 

n— 1 

the number of trees in A„ is # m (n) = ((to — 1 )k + 1). 

k =1 
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3 The dynamics 

Let N be the (large) number of interacting components. Let m (m > 2) be the 
fixed number of components involved in each interaction. We suppose that all 
components take their values in a measurable space, (E,£), (one can think of 
(R d , B(M. d ) or simply a finite set) and their interactions are given by a symmetric 
probability kernel Q on the product space (E m ,£^ m ) . That is, the function 
Q(x\, X 2 , h-Xm'i C i x • • • x C m ): is measurable in (xi,X 2 , ■■■x m ); is a probability 
measure in {C i x • • • x C m ); and satisfies Q(x\, X 2 , ■■■x m ',C\ x • • • x C m ) = 
Q(®<7(1),®<7(2),—*<T(m);C'<T(i) X C<t(2) X • • • x C^ m )) for any permutation a of 
{1, 2,..., m}. 

In the following example, we simplify the model of Duffie-Garleanu-Pedersen 
[6] by keeping only their binary interacting kernel. 

Example 1 . Investors in this model have two liquidity states denoted h, for 
high, and l, for low. Moreover,there is an asset of common interest to these 
investors who either own the asset (denoted by o) or don’t (denoted by n). So 
E = 0)} describes the state space. The kernel is defined 

by Q 2 )-, •; C\ x C 2 ) = 0 except for 

Q 2 ((h, n), (l, o); Ci x C 2 ) = Q 2 HI, o), (h, n); C 2 x Ci) = 5(h,o){Ci)\i, n ){C2) 

where S Zo is the Dirac function S Zo (z) = 1 iff z = Zg and S Zo (z) = 0 otherwise. 
The binary kernel implements the trading of the asset whenever a low liquidity 
investor who owns the asset meets a high liquidity investor who does not yet 
hold it. 


The interactions occur at each jump of a Poisson process with intensity A^. 

Groups are undistinguishable so each group has a probability of (^) 1 of being 
involved in a given interaction. 

The kernel Q allows us to describe the macroscopic evolution of the system 
with an associated system of non-linear differential equations via the evolution 
of the law of a component. This probability law, denoted p t , evolves with time 
and is in fact the solution of the Cauchy problem: 


dut 

dt. 


= Hl^t 


Ht) ;Ak) = h 


where 


p° m (C) = jp(dx\)p,(dx 2 )-..pL{dx m )Q(xi,X 2 ,...x rn \C x E m 1 ) for C 


G £. 


The probability law y,° m is the law of a component after the interaction of 
m i.i.d. components with law p. We can think of it as the law at the root of the 
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first meeting 


time ti 

second meeting 


time t? 


Figure 1: Simple interaction tree with only two meetings and m = 3. 


m -ary tree with only one interaction. We will look at all the trees representing 
the interaction history of a component up to time t. So for a tree, A, with more 
than one interaction, we divide the tree in m subtrees at that last interaction 
and continue recursively up to time 0 to define \x° mA . (Please see figure 1 for 
a simple example of an interaction tree.) Let A„ be the set of all trees with 
n interactions (a.k.a. nodes), each node producing m branches. If A n £ A ra , 
then denotes the law obtained by iteration of fj,° m through the successive 

nodes of the tree when we place the law /i on each leaf of A n . 


We have shown Belanger-Giroux [1] that the Cauchy problem has a unique 
solution which can be expressed, by conditioning on the number of interactions 
up to time t, and then by the component’s history. Such conditionings give us 

Vt = ^2Pn(t) n , s E ^ W 

«>0 fra " j A„£A„ 
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n— 1 

where # TO (n) = n ((to— 1)/c + 1) is the number of trees with n nodes, taking 

fc=1 

into account their branching orders; and p n {t) = e~ A *(l — e -( m - 1 ) At ) n 

is the probability of having n branchings up to time t. 

Remark 2 We call the law p t = y^e~ At (l — frri _] )rln] ^ p ° mAn 

n>o A n eAn 

an explicit extended Wild sum [15] and note that the convex combination we ob¬ 
tain for the case m = 2 is indeed the Wild sum, pt = ^^e _At (l—e _At )"^r p° mA 

n >o A n e A„ 

now well-known in the statistical physics of gases since the work of Kac (1956) 

[ 8 ]. 

3.1 Using interaction trees to go from the microscopic to 
the macroscopic. 

In all our cases, we have an underlying market structure which is a Kac walk 
with interactions involving to agents. We add exponential times to obtain a 
marked Poisson process whose marks are horizontal lines linking the agents 
participating in a given interaction. This enabled us, in Belanger-Giroux [1], 
to describe the limit law of an agent, under an appropriate conditioning, as a 
countable convex combination on trees which is, as we have shown in section 3 
of that article, the global solution of the associated differential equation on the 
space of probability laws. 

Here we first explain how we came to that convex combination since it serves 
as a tool to study the other models which follow. It is the tool that enables us, 
for instance, to state proposition 5. Its proof follows the lines of the proof of 
the main result in Belanger-Giroux [1], 

We start our study by an analysis of the dynamics of the intrinsic structure 
of the large set of interacting agents when the number of agents increases. We 
assume that each interaction involves m agents, in > 2. More specifically, 
we consider a set of N agents whose interactions happen at unexpected times 
so these interactions’ occurrences follow a Poisson process. Since agents are 

interchangeable, each group has an equal probability of meeting of 

If we suppose the intensity of the meetings to be ^ then each agent has a 
meeting rate A which can be assumed to equal 1 under a time change. Wc will 
make this assumption, A = 1, all throughout section 3. 

For N fixed and starting at time 0, we assign a vertical position to each 
agent. The down movement represents the passage of time, see figure 1 on page 
4. Each time a group of agents interacts, we draw a horizontal line between 
those agents and we draw a vertical line at each agent’s position connecting 0 to 
the horizontal line just drawn, so we see a random graph being formed. When 
we stop this graph at time t, we obtain the finite graph of all interactions that 
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first meeting 


second meeting 


time t 
P 

Figure 2: Simple graph with cycle. 


have taken place. Moreover, the history up to time t of a given agent, call it 
P, is described by the random graph connecting all agents who have interacted 
directly or indirectly with P. 

The number of meetings is random but we can condition on it. The law of 
the finite graph is reversible since the meeting times are uniform on [0, t ]. We 
want to show that a random graph representing the history of P can be replaced 
by a random tree as the number of agents, N, grows. If we look at figure 2, 
we see that the inclusion in the second meeting of one of the investors having 
participated in the first one would create a cycle in our graph. As N grows 
though, the chance of meeting an investor previously encountered directly or 
indirectly tends to zero. 

To see this, let us consider the graph of P’s history up to time t. Starting 
at time t, we pursue each one of the encountered vertical lines in P’s history 
backward in time until we reach the next horizontal line. If the inclusion of the 
horizontal line in our graph does not create a cycle (i.e. no pair of investors 
were involved directly or indirectly in a previous meeting) we include the line, 
if not we remove it. Proceeding in this fashion up to time 0 we get a tree with 



n internal nodes, say, which has the same law as the law of a tree obtained by 
a pure-birth process. The tree obtained by a sample history of P’s interactions 
is an ?n-ary tree. 

These trees grow randomly in time: each time a new node appears, corre¬ 
sponding to the occurrence of a meeting of investors at that time. We recall 
that A„ denotes the set of ro-ary ordered trees with n internal nodes. Then A„ 
constitutes a set of random trees if we assume that every to- ary tree in A„ is 
equally likely, namely of probability # 1 ^ n ~ ) . 

The tree starting at P's vertical line at time t with intensity 1 and which at 
time 0 has intensity (to — l)n + 1 and that same number of leaves. Between two 
branchings of this process a graph representing P’s meeting history can have a 
random number of additional horizontal lines following a Poisson law of param- 
(to — 1 )n + 1 \ / N 
2 M TO 


eter at most ^ 


-IN 


We will now bound the ex¬ 


pectation of these supplementary horizontal lines by a majorant which tends to 
0 as TV increases. Indeed, since the mean number of redundant lines when there 

are n branchings up to time t is at most I (( (m “ ^ + 1 ) ( ^ 

we have that the mean number of redundant horizontal lines is bounded above 
by 

/ (?n — 1 )n + 1 \ ( N 
m \ l 2 J \ rn 



PN,n(t ), 


where PN,n(t) is the probability of having n branchings up to time t of the pure 
birth process with successive branching waiting times following exponential laws 
of parameter 


N 

A N,n = — ((to - 1 )n + 1) 
TO- 


N — ((to. — l)n + 1) 

TO. — 1 


N 

TO. 


Since 


N 

— ((to — l)n + 1) 


N — ((to — l)n + 1) 


TO " ' \ TO - 1 

_ ((to. - l)n + l)( Ar - (( ™I 1 1 )n+1) ) 

C=D 

< (to. — l)n + 1 


N 

TO 


-1 


( 2 ) 


then pN,n(t) is stochastically smaller than the law obtained with the intensities 
X n = (to. — l)n + 1, which in turn are less than the intensities X n = m(n + 1). 
Its transition kernel is then obtained by solving Kolmogorov’s affine system of 
equations: 
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mp t { 0) 


dp t { 0) 
dt 

dpAn ) _ . . . . 

——— = mnp t (n — 1) — m(n + 1 )p t (n) ; n > 1. 

Thus the latter intensities give us a geometric law p t {n ) = e~ mt (l — e~ mt ) n = 
e -m(n+i)t ( e mt _ _ gi nce geometric laws have finite moments of all orders, 

the mean number of redundant horizontal lines is bounded above by a quantity 
converging to 0. 

For more details on Kolmogorov systems of equations for pure birth processes 
we refer the reader to Lefebvre [10], for instance. 

Thus, after having specified the initial agents’ states and their interaction 
kernels, we can approximate P's law using the tree obtained from removing all 
redundant horizontal lines from its graph. We will use this fact in the next 
sub-section. 

3.2 Limit countable convex combination 

We will now show that these random trees whose branching intensities depend 
on N can be approximated by trees with branching intensities independent of N. 
Taking into account that P’s tree history is random with intensities depending 
on N, we could write P’s law, denoted by pf :N , with complex formulae depend¬ 
ing on N. Since our markets have a large number of investors, it is preferable 
instead to work with the limit of these laws. We note from (2) above that for 
each n, Xn.u —!► ((m — l)n + 1) as an increasing sequence in N. 

Let p n {t) (— Pt{n)) be the solution of the affine Kolmogorov system of equa¬ 
tions: 

^-,(0) (3, 

= ((m - l)(n - 1) + 1 )pt(n - 1) - ((m - 1 )n + 1 )pt(n) ; n > 1. 
dt 

Recall fron the first section that p t = Z! p ° mAn - 

n>0 A„ 

Proposition 3 The sequence of laws converges to ji t as N increases. 


Proof. By Kurtz [9], we have that PN,n(t) Pn(t) as N increases. But 
(Pn(t)) n >o is a probability law, so for e > 0, there exists n(e) such that 

Pn(t) < e. Now let N(e) be such that N > N(e) implies that \pN,n(t) — 

n>n(e) 

p n (t) | < f° r 0 < n < n(e).We then have for C G £ and N > N(e) 
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n(e) 

|/V W (C) - Ht(C)\ < ^2\p N , n (t) -Pn{t)\ +2e < 3e 

n =0 

since #^Xn) M 0 mA "(C) < 1 and ( PN,nit )) n >o are probability laws. Our 

.4„eA„ 

claim is proved. ■ 

Lemma 4 p t (n) = - e^™- 1 ) 4 )” 


Proof. We need to solve the affine Kolmogorov system of equations (3). 
Proceeding by induction we have: 


dpt{ 0 ) __ t 

dt 

t 

^4^- = ((m - 1 ){n - 1 ) + 1 ) e -(n(m-l)+l)t r e ( n(ro -l)+l ) 8 ( n _ 1)da 

dt J 

o 

To prove the lemma it suffices to note that # m (n) = # m (n—l)((n—1 )(to— 1) + 1) 
and that e (”( ro - 1 )+ 1 ) s e“ s (l - e -(m-i) S )n-i = e (m-i) s ^ e (m-i) s __ ^n -1 is the 

derivative of {m \ )n (e (m ~ 1)s - l) n . ■ 

And this shows that the limit law of P is indeed the extended Wild sum 
which we have shown (in [1]) to be the solution of the ODE associated to the 
interacting system. 


4 Explicit formulas for other OTC market mod¬ 
els 

In many applications, it is more convenient to work with more than one kernel 
to describe the dynamics of the system. It is the case for instance in the models 
of Duffie-Garleanu-Pedersen [ 6 ] and their extensions in Belanger-Giroux-Moisan 
[2] and in Belanger-Giroux-Ndoune [3]. 

In the simplest such model on E = 0)} we have the 

binary kernel we described at the beginning of section 3 and we have the au¬ 
tonomous changes of liquidity of an investor. Let j u and "fd resp. be the intensity 
of the up movements (resp. down movements) in liquidity. We will first assume 
that these intensities are equal (we will remove this assumption at the end of 
the section) and we let 7 = = 7 d- Then q p (t) = e _ 74 ^j- is the probability 


9 



of having p autonomous movements up to time t. The 1-ary kernel can then be 
defined by 

Qi((l,n)-,C) = S( hin )(C);Qi((l, o);C) = 6 {h , o) (C); 

Qi{(h, n);C) = 5 M (cy,Q 1 ((h,oy,c)=s m (cy, 

and 

v 01 (C) = i/(l, n)S( ht „) (C ) + i/(l, o)6( ht0 ) (C) + i/(h, n)S^ >n ) (C) + u(h, o)S ( i t0) (C) 
for C C E. 

It is possible to modify the 1-ary kernel into a binary kernel with the addition 
of a “witness” investor who is completely unaffected by the change of liquidity of 
the other investor. We then symmetrize that kernel and replace the two binary 
kernels with a convex combination of the two kernels to be left with only one 
kernel as in the situation we dealt with in the preceding sections. So all the 
formalism developed so far is still valid. The drawbacks of this approach though 
are that the synnnetrization operation gives us a slightly different dynamics, and 
more importantly, that we lose the explicit formulas for the extended Wild sums. 
The objective of this section is to show how we can retain them. 

Let K p denote the set of all arrangements of p undistinguishable objects in n 
boxes, where a box may contain arbitrarily many objects. Then |K P | = ( n+ ^7 )■ 
Let A n denote, as before, the set of all random trees with n m— ary interactions 
(the investors meetings). If A n £ A„ then A n has (m — 1 )n + 1 leaves which in 
turn gives mn + 1 branches. 

Let p denote the number of 1—ary interactions (i.e. the number of au¬ 
tonomous changes of position). For a G K^ n+1 , let A a n denote the tree obtained 
by placing the 1-ary interactions on each branch of the tree according to the 
arrangement a. Let A rlJJ denote the set of trees with n in- ary interactions and 
p 1-ary ones. Then p : K^ n+1 xA„-> A„ iP : (a, A n ) i-G A° defines a bijection. 
Please see figure 2 for simple examples of trees in A 2 ,i- Moreover, if we call 
{<7*}J =1 the 7 configurations of figure 3, then (i/° 3 A2 )°i = ^J =1 v ° sA 2 . 

The set A„ iP of (m, l)-ary ordered trees with n m-ary internal nodes and p 
1-ary nodes has the cardinality equal to # m (n) ( m ”^ p ). Then A„ iP constitutes 
a set random trees if we assume that every (m, l)-ary tree in A„ iP is equally 
likely, namely with probability — —^y^rs+pj- 
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P(Aj)=/ t_ 


Figure 3: Interaction trees for m=3, A: and p=l. 


Proof. 

Let /It. be the solution (1) above. Since \vt\ < |/it| and fi t is uniformly 
summable, the convex sum v t can be differentiated term by term to obtain: 


dvt 

dt 


-A v t + 


Xe 


—mXt 


Ea-< 


— (m—l)Xt\n—l 


1 


n> 1 


(m — l) ra_1 (?r — 1)! 


E 

A„eA„ 


E gpft) o m A 

/ mn-\-p\ / j 

p= 0 V mn ) crgK p , . 


-T^ + E^W^n) ^ 

n >0 4 




P - 1 


v c -‘(7 ty 


E 


,°mA" 


■6K P 


Now using the combinatorial results of the proof of theorem 1 in Belanger- 
Giroux [1], we have that vl™ is equal to the expression 


Ae -mAr^(i_ ( 


— (m—l)Xt\n—l 


a>l 


(m — l) n_1 (n — 1) 


1 (n — 1)! .E 


A n e A„ 




^ ] /mn+p\ 
p =0 V mn / crgK p 


E - 


,o rn A 


Otherwise, the expression 

Ep»( % x , \ E 

n>0 ttm( n ) Ane A „ 


v e _ v ^ 

iq/mn+m ^ 

P=1 V ^ ''•V mn 


o m A p 


is equal to the quantity 


n>0 


1 


E 

:nGA 


v- e 7 (7^) p Y' i, 0 ™- 4 ™ 

Z^„!('mn+p+ 1 \ 

P=0^A mn ) ffgK P+i +i 


But ^ 

aer+' + i 


= E ( v° mA "V 1 .And then 

<r'eIK p ,, 
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^2Pn(t) „ 1 , , ^ 

n>0 ttrn{n) An&Ari 


y e v (j/ 

Z^j/mn+p+n Z-f ' 

P=0^'V mn / cr'glRP 


,o m A" \o 


) 01 


Hence we get that ^ has the desired form and the proof of the proposition 
is now complete. 


We note that if p = 0, that is, if no investor changes its liquidity position, 
the above solution does indeed become the solution (1). 

Remark 6 In the specific context of the DGP model we have a binary kernel 
which simplifies the first part of the formula. But without the assumption j u = 
"fd> we have to consider the up movements and the down movements separately, 
and this makes for a more complicated second part of the formula. Let 7 = 
7 m + 7 d, the solution becomes 


v t 


= E- 

n> 0 


e (l — e _t ) 


E 

A„eA„ 




o-7* 


mn+p+l\ 


E(^) fc (7 rf ) p - fc E 


o A n 


k—0 


^ slK m _ n+1 


where a u (resp <7d) denotes the arrangements of up movements (resp. down 
movements) on the branches of the tree and a u Uad is the arrangement obtained 
from both arrangements of up and down movements. 


Remark 7 We can obtain similar explicit formulas for OTC models where the 
interactions involve m > 2 investors. In the information percolation model of 
Duffie-Malamud-Manso[7], for instance, the state space, E = N represents the 
potential levels of information acquired by an investor through meetings with 
other investors. The m-ary interaction is the perfect sharing of information 
which means that each investor in the meeting comes out with the sum of the 
information levels of all participating investors. The unary kernel is a regression 
force which replaces an investor of level n say, by an investor with level n(n) 
sampled from a given distribution tt on N. 
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